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Abstract 

Motivated by the similar mass splitting in light-light and heavy-light = 0~ and = 0^ mesons, 
the SU(3)-symmetry breaking effects splitting the massses in the 0~ and 0"^ channels of the D meson 
are analyzed in the framework of QCD sum rules with an underlying cq structure. We take into 
account operator mixing to obtain an infrared stable OPE including complete non-perturbative and 
perturbative 0{mq) corrections to the correlation function. With the same threshold for both channels, 
the mass splitting arising from the sum-rules has the same behavior as the observed spectrum. In 
particular, we obtain mOs — ~ 35MeV in the 0~ channel and mo^ — mos ~ 12MeV in the O"*" 
channel at a renormalization scale /.t = IGeV. The splitting can be attributed to the different roles of 
mass effects and the parity-dependent " force " induced from non-perturbative QCD vacuum. Further 
analysis shows that due to this "parity-dependent" force it is natural that the mass gap of the two 
states in the 0" channel is larger than the 0"*" channel. When we increase the renormalization scale to 
IX = 1.3GeV the splitting remains unchanged which demonstrates a correct scale invariance. Combined 
with HQET, generalization to other channels of charmed mesons and 6-systems are briefly discussed . 

PACS numbers: 12.38.Lg, 14.40.Lb 
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INTRODUCTION 



m 



The SU(3) quark model of hadrons[l|, |2| provide an intuitive understanding of hadronic 
properties. Due to the non-perturbative nature of low energy QCD we have to employ non- 
perturbative methods in the hadronic sector. The QCD sum rule approachjs-Sl has proven to 
be a successful non-perturbative method to extract reasonable results in the hadronic sector. 
Similar success have been achieved with light-cone QCD sum rules which represent a 
further-developed version of the original sum-rule approach. We will not dwell on the overall 
success of the QCD sum rules, but will focus on the light and heavy = 0~ and = 0"*" 
channels. 

361 ] that the mass splitting in the = 0~ channel for the lowest light mesons is 



It is observed 



in line with their underlying structures from the naive quark model estimates. This splitting of 
lowest light pseudoscalars with quantum numbers of the vr, and rj' is well accommodated 
in QCD sum rules if the instanton effects are appropriately |8| included since the instanton 
contributions to the correlation function are different from each member of a multiplet due to 
its dependence on the isospin and effective mass m*^. However, the splitting in light scalar 
meson = O"*" is the reverse of the naive quark model estimate. If instanton effects are 
considered, the splitting in the = 0+ channel above IGeV (i.e. /o(1370), ao(1450), A'o*(1430) 
and /o(1500))can also be explained within the framework of the QCD sum rule approach^. 

If we assume an ordinary light-heavy underlying structure of open-charm systems a similar 
mass hierarchy as light pseudoscalars can also be observed in the = 0~ channel of D mesons, 
the L'd(1869) and 0^(1968). HoweverMn the = 0+ channel L),(2317l (which was first 



discovered by BARBAR Collaboration 



10| and later confirmed by CLEO ll|]) and its isospin 



partner D*{2A00)~^ (observed by FOCUS Collaboration 12|) also show similar splitting as light 
scalars of = 0"*" channel in contradiction to the naive quark model estimate. Among these 
open-charm systems Ds{2317) triggers much attention on its underlying structure. Mass results 
from Lattice QCD for Z}s(2317) are larger than the experimental value 13l-ll5| and the results in 



Ref. 



13| suggested that Ds{2317) might receive a large DK component. The work of Ref. 16] 



including this contribution from DK continuum in QCD sum rules based on a cs structure 



^ One should not confuse it with the " effective mass " in the mnsituent quark model. 
^ The results show there is large glueball content in /o(1500).[9[ 



2 



found that this continuum contribution can significantly lower the mass and decay constant of 



the Z^s(2317). A more complete work on open-charm systems can be found in 17| where the 
JP = 0-,0+, channels were studied from a viewpoint of cq system where q = u,d, s. 

On the contrary, in considering the difficulty of ordinary heavy-light structure in decoding 



the nature of Ds{2317), a four-quark states picture was proposed[l9|, l20|.Ref. \2l\ employed 
the four-quark structure to investigate Ds{2317) using QCD sum rules which suggested that 
Ds{2317) might be a four-quark states, while the radiative decay of i5s(2317) in light-cone sum 



rules favors a cs structure (22l|. All the work both experimentally and theoretically shed some 
light on the interpretation of Ds{2317). 

One might hope that the instanton improved QCD sum rules can also realize the splitting 
in pseudoscalar and scalar D mesons. But this seems unlikely because there will be rapid sup- 
pression of instanton contributions by the c-quark effective mass m* and damping exponential 
factor in the case of the charm quark, and thus the QCD vacuum condensates of various opera- 
tors are still the dominant non-perturbative corrections. These condensates do not preserve an 
ideal SU(3) flavor symmetry because of the symmetry-breaking effect of different quark masses. 

In the framework of QCD sum rules the SU(3)-symmetry breaking effects present itself 
mainly from two sides. One is the different values of (uu) , (dd) and (ss) , and associated mixed 
condensates. Ref.j23| studied the role of the chiral condensate (qq) in the mass splitting be- 
tween the scalar-pseudoscalar D(and B) mesons with so-called no-free parameters sum rules 
in the chiral limit, and recently a parallel analysis was applied to the dependence of heavy- 



baryons mass splitting on the ratio k = {ss) / ((id) j2J]. Another symmetry-breaking effect is the 
perturbative mass correction which is proportional to mq^mq^ where qi and q2 labels the quark 
content of the meson or the current considered. These breaking effects are small for pure-light 
mesons since the quark condensates, mixing condensates as well as the perturbative mass cor- 
rections are always accompanied by the light quark masses and therefore they will be greatly 
suppressed, especially for and m^. In fact we always use a massless approximation for pure 
light systems. However, when heavy quarks are involved, the mass effects will be considerable 
since the large mass of heavy quark takes the place of one of the light ones mentioned above. In 
other words the large mass of the heavy quark in heavy-light systems results in more signiflcant 
mass effects than the pure-light system. Thus it is more consistent to take into account the 

n 

mass corrections to an uniform order. The workp/7| went further by including the perturbative 



mass corrections, and the splitting in different states was realized by choosing different thresh- 
olds. However, the theoretical splitting(given by the central value) in 0"^ channel, although 



in agreement with the potential model results 18|, still contradicts experiment. Therefore a 



complete analysis including both operator mixing effects and all the dimension-6 operators is 
necessary. 

Motivated by the important role of heavy quark, in this work we will use QCD sum-rules to 
investigate the mass effects in splitting the pseudoscalar and scalar D multiplet. To be specific, 
in our work complete perturbative and no-prerturbative 0{mq) mass corrections (where 
is the mass of light quark) are taken into account in the sum rules. Unlike pure-light meson 
SU(3)-breaking effects mainly introduced by VEVs of renormalization-group invariant mass- 
dependent operators, in heavy-light system there will be another source introduced by operator 
mixing mentioned above. Another important point is that in QCD sum rules the threshold also 
represents SU(3)-symmetry breaking, so different members of the same SU(3) multiplet should 
have different thresholds. In practice we have no knowledge about SU(3)-symmetry breaking 
effects a priori; it is only evident from experiments. For instance we even do not know how 
to set the thresholds for the 0"'' doublet before we know their exact masses experimentally. 
Normally if we take the viewpoint of constituent quark, the threshold of Ds should be larger 
than that of Dd- But if the mass of is larger than that of Dg in O"*" channel it would be not 
reasonable to set a larger threshold for Dg. For this reason we would like to focus our attention 
on the splitting trend rather than the exact spectrum of the D mesons. We believe the mass 
splitting of the doublet Dg and D^, in QCD sum rules under the same "reference" threshold 
reflects the SU(3)-splitting tendency. Thus for an exact sum rule a suitable modification of 
threshold due to SU(3)-breaking should enlarge this tendency, otherwise we think it is not 
natural. To this end in our analysis we first select one member of a multiplet as " benchmark " 
to fix the suitable threshold and Borel window. Then we apply these parameters to another 
member with replacing the SU(3)-breaking dependent quantities. 

The article is structured as follows. In Section II we first review the necessary results on 
the operator mixing and cancelation of mass singularities for heavy-light current, then present 
the sum rules for pseudoscalar and scalar currents of mesons. In Section III the numerical 
results and discussion will be given. Conclusions are presented in Section IV. 
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II. THE FORMULAS 



A. operator mixing and cancelation of mass singularities 

In order to demonstrate the operator mixing and cancelation of mass singularities in heavy- 
light quark system, we consider the following charmed scalar two-point function: 



n(g2) = i j d4xe*9x^o|T{g(a;)c(x),c(0)g(0)}|0) 

= C,(g')/ + $^C,(g2)(0|O.|0), 

where d is the dimension of the operator and q = u,d, s. For simplicity the renormalization 
invariant factor (ln(yu/A))~^/'' has been suppressed where /i is the normalization point and 
b = {llNc — 2nf)/3. Setting aside the perturbative part until later, the contributions of VEVs 



of d < 6 may be written as 271] 



71 71 



[ 2 71 12m. 

+ {c^q) 

24 /i2 



+Q,,|(gjg)-^ln^(j^)^ (1) 
The operators in Eq.(IT]) are defined as follows: 



71 \ 7T 



I - \9sJ ^,„y^,yp^p^. 



if) = (^g',{D,G^^niD^G;^)) = U{Y,qYT'^q 
{qjq) = (gsqi^iD^^nT^q) = lgl[qi,T-q) J^q^T'^q 
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therefore the contribution from (j^) is to 0{al), in the following we assume the singularities in 
are canceled by the mixing and we will omit this term. All the charmed-condensates vanish 
by virtue of heavy quark expansion [sl. I271 29 1 : 

1 ,«« „o, 1 



(cc) 



(cGc) = ^ In ^(^G^)- 



12mc TT 
2 



1 ,a 



G') + ... 



-G-') + ... 



(2) 

thus there are only gluonic and light quark related condensates left. It is clear there are mixing 
to gluonic operators from (qq) and (qGq) terms, and with the help of this mixing the singular 
parts in gluonic coefficients in limit nig — will be well canceled and we are left with an infrared 
stable expression. The final- form of various non-perturbative coefficients follows Eq.([T]): 



Cgq 


— r 


Gg^ 


= Cg^ + 


GqGq 


= CqGq, 




= Cqjq, 


Gg3 


= Cgs + 



U^q^'' - ^ ill 



360m3 Uruq 



G 



qGq- 



Three of the Gs have been worked out 



30| in expansion in niq. In our notation: 



G 



rur. 



m 



g2 _ ^2 



+ 



2 s 



TUq 2m^ — q^ 
2 (g^ — m^)^ (g^ — m^) 



2N3' 



12(g^ — ml 



nir 



niq 2(g^ 



+ 



+6m^ ( - In — + In ^ 



1 , , /imc 



C. 



^2 1 ■'"<? 
.2 /i 

,2 



qGq 



niqml 



mc{q^ — 



"im m^^q" m 
In — - 



q + 6m 
2 



2 (g2-m2)3 



2 (g2-m2)4' 
and the singular pieces of Cg'^ as m„ — )■ are^ p 



Gg3 



nrir 



+ 



g^ - 2ml 



360mj^(g2 — m^) 720m2(g2 — m^)^ 



mc(15g^ — m^) 
360mq(g2 — 



2^3 



(3) 



(4) 



(5) 



^ The three-gluonic coefficient in 28|was derived from pscudoscalar heavy-light current, here one can obtain 
the scalar one by replacing rUc by —mc- 
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Substituting Eq.(|l]) and Eq.(|5]) into Eq.Q all the mass singular parts of Cc-a and C^s appearing 
as l/mq and Inm^ are canceled since these terms are remnants of long distance structure of 
vacuum condensates. For definiteness we write down the explicit form of (7^2 to 0{mq) as: 

1 mr,m^ 



Co2 



3 



12(g2 _ m2) 12(g2 - m^) 

+ !M! ^ ("gS ^ Q^2 ^ g^2 _Jfn^\ ^ /gv 

12mc(g2 — m2)3 V 77^2 _g2y' 

Therefore to 0{mq) the operator mixing changes the C^a significantly and it is expected there 
will be new mass effects on the sum rules. Similarly, for Cg-s we have(see the appendix for 
details) : 

Cg^ = -7^^(-10^' + + + 2), (7) 

where 

W= 



ml — 



Now we have fixed the non-perturbative parts in the OPE of scalar heavy-light current. It is 
easily to get the non-perturbative parts of pseudoscalar current by replacing mc by —m^ in 
Eq.Q. The perturbative part and Cqjq will be presented in the forthcoming subsection. 



B. the sum rules 

The sum rules of scalar and pseudoscalar D are based on the following two-point correlation 
function: 

n.„=)^./.Se-(0|r{«wrcW,.,0)r,(0)}|0>. (8) 

where q is the light flavor in the D meson, F = {1,275} for scalar and pseudoscalar D meson 
respectively. The decay matrix element of scalar D meson is defined as: 



(0|gc|0) = mz)/. 



D-, 



and following 33| the pseudoscalar one is defined as: 

mc{Q\qii^c\D) = ml)f'j^, 

where m^ is the c-quark mass and m^, the D mass, and labels the mass of light quark. In- 
cluding the perturbative mass correction to 0{mq) and to 0{as) as well as the complete 0{mq) 



non-perturbative corrections worked out in previous subsection, after Borel transformation the 
OPE of correlation function in Eq.dH]) is given by^ 3ll-l33l|: 



no™ 



2 

mi \ 




f{s,m 



exp[ 



M2- 



\i ,1 ^"^^ 



4M6 
27M4 



4M2 12M4 



(9) 



where 



with 



and 



f{s,m 



9 
4 

+ ln 



2Li2 



In — - In 



3 , m 
-In : 

2 s 



3 ml ^ s — ml 

— ^ + ^ In ^ + - 

mi s mi s 



m 



2 ^'^—^^ 
mi mi 



Li2(x) 



dt 



ln(l -t) 



7^ = 0.577. 

is the Euler constant. Contributions from 3-gluonic condensates have been omitted safely 
since it is m^-independent and greatly suppressed by the huge denominator therefore it is not 
responsible for the splitting as one can see from Eq.(I7]). The upper and lower signs in Eq.(|9]) 
are for the scalar and pseudoscalar channel respectively. 



The mg-independent term in coefficient of (^G^) here is different fromflTt and 23 1, but this term has no 
impact on sphtting except an uniform shift. 
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On the other hand the correlation function in Eq.([H]) can also be derived from the phe- 
nomenological side by the dispersion relation: 



as h subtraction constant, 



(10) 



where the spectral density Imnp^(s) is obtained by inserting a complete set of quantum states 
S|n)(n| into Eq.® which reads: 



Imnr(s) = FrTn^j^n 5{s — mjj) + n- 



3 (s — m 



2\2 



It 



2mqmc 4as 



s — mi 



+ -±f{s,mi) 6(3 -So), (11) 
on 



Taking the Borel transformation of Eq. (fTO|) and equating it with Eq.([9]), after subtracting the 
continuum contributions we arrive the desired sum rules: 



Frm'jj exp[ 



m 



87r2 



ds — 



+ 



12 ^ 



M2 



s — rat 



mq{qq) 



exp[ 



exp[ 



12mr 



8M4 



2mi 



2M2 



M2 



M2 



(^G^.G-)exp[-^] 



1 / 777^ \ \^ 

± ^( 1 ^-)mJqga.^—G^^'^q) 



4M4 

lOvra^ 
27M2 



4M2 12M4 



{qq) }exp[ 



M2J 



(12) 



Now we have completed the sum rules for pseudoscalar and scalar D mesons. The input 
parameters in Eq. (fT2|) are as follows 



34 



36 



4Q|: 



Aqcd = 259MeV, a, = 0.517, (a.G^) = (0.07 ± 0.01)GeV^ 
{uu) = (dd) = -(0.225 ± 0.025)^GeV^ {ss) = 0.8(nn), 



mrf = 3 ~ 7MeV, 



m, = 120MeV, 



2/^ 



2 = O.SGeV^. 



(13) 
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All the values adopted above are given at the scale u = IGeV and we deduce the QCD scale 



Aqcd to one- 



is given by|35|: 



oop from as{Mz) = 0.1170 ± 0.0012[26|. The renormalization scale dependence 



a. 



TT 



TT 



(14) 



with b = (llA^c — 2?7,j)/3. We use the following pole mass for the charm quark^ 2q : 



= 1.47 ± 0.04GeV, 



which can be expressed in terms of the running mass through the relation: 



mc = mc(/i) 



'+ln4l^ + 0(«^) 



(15) 



^ 3 / TT 

Taking the logarithm of both sides of Eq.( |T2l) and applying the differential operator 
M^d/dM'^ to them we can separate the mass from decay constant. Now we have fixed all 
the ingredients for numerical analysis. 



III. RESULTS AND DISCUSSION 

Firstly we present the criteria followed in our analysis: 

1. To specify the appropriate threshold and Borel window, we demand that the continuum 
contribution [i.e., the part in the dispersive integral from sq to oo which has been subtracted 
from both sides of Eq. f|T2l) ] should not be too large (less than 30% of the total dispersive 
integral). This criterion give us an upper limit on the Borel momentum M^. Furthermore, the 
non-perturbative dimension-six operators corrections should be less than 10% which establishes 
a lower limit. 



^ A note in 23| argued that the value m,c = 1.46GeV used in[17[ may be iU-defined. As the sum rule is sensitive 
to nic as we can see in the next section, the larger choice rric = 1-47 ± 0.04GeV here might induce new error. 
But it does not affect the splitting since the thresholds are fixed for same channel in our analysis. 



10 



2. In order to check the mass effects on the sphtting of the same channel, states of that 
channel will be analyzed under same threshold while the flavor-dependent parameters such as 
TTig, {qq) as well as (qcrGq) change correspondingly. This will supply us with an appropriate 
comparison in the same channel with different light content. 

3. As mentioned above our primary concern is a correct splitting trend, not the whole 
spectrum in same channel. Thus we select one state as our "benchmark" to determine the 
threshold and Borel window according to criterion 1. For definiteness we select and Dg as 
our sample in 0~ and in 0^ channel respectively. After this we turn on another state following 
criterion 2. If it is a natural sum rule, it should produce a correct splitting trend that agrees 
with the experiment. 

With these criteria in mind we plot the mass curves of the two states in the same channel 
against the Borel momentum in a diagram for different threshold and charm mass rric since 
it is convenient to observe the splitting. The working windows which satisfy the criterion 1 
are marked by two short lines(or one short line which labels the upper limit only) while the 
narrow ranges from which we read our numerical value are marked by shaded bars. If there is 
no an obvious extremum within the window we determined, the central value will be adopted. 
Under these criteria we find for fixed threshold and charm mass rric as well as scale parameter 
fi the working windows for Dd and Ds in each channel are very close. The upper limit of 
our working windows decrease as the thresholds decrease, while it seems that the lower limit is 
nearly invariant which can be seen from the following graphs. When we scale up to /i = 1.3GeV 
the upper limit increases compared with fi = IGeV while there is no obvious impact on the 
lower limit. 

First, we study the 0~ channel from the pseudoscalar sum rule given in Eq.(!T2l) at the scale 
/i = IGeV. The numerical results of the two states of = 0~ channel are shown in figureJH We 
can see from figurelUthat the results following the pseudoscalar sum rules of Eq.f ll2p accurately 
reproduce the mass of with a large charm mass which is very close to the experimental 
Z)rf(1869). When the parameters of Dg turn on and with the threshold fixed, the resulting mass 
from the sum rule Eq. f|T2|) is still lower compared with the experimental one Dso(1968) but it 
is always larger than Da- The mass splitting from Eq.( !T2l) is ~ 35MeV, a value much lower 
than the observed splitting ~ lOOMeV. At a first glance it appears that the SU(3)-breaking 
effects can not supply a realistic splitting. But it is important that the sum rule does present a 
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2.08 - 




0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2.0 

M''(GeV') 



FIG. 1: Mass of ^'^(dashed line) and Dg meson (solid line) of channel from pseudoscalar sum rule 
of Eq. (jl2p as function of Borel momentum Af^ at scale ^ = IGeV. 

reasonable mass splitting trend in the 0~ channel which is in line with the observed spectrum. 
The failure of pseudoscalar sum rules of Eq. f|T2|) in providing the entire mass gap in 0~ channel 
is understandable: there is a large gap between the two 0" states thus a threshold appropriate 
to 1^(^(1869) is too low to produce i5s(1986). On the other hand in QCD sum rules we notice 
the difference of threshold between different members also reflects SU(3)-symmetry breaking. 
Therefore it is expected that if we determine the threshold of each member separately, at some 
lager threshold than the one for Dd the mass of Dg will be well produced from the sum rules. In 
addition, the theoretical results are very sensitive to charm mass: as we can see from figureiH 
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1.7 1.8 1.9 2.0 2.1 



M'(GeV') 



FIG. 2: Mass of (dashed line) and Dg meson (solid line) from scalar sum rule of Eq. ()12p as function 
of Borel momentum at ;U = IGeV. 

at the lowest value adopted in our work, even at a much higher sq it is still difficult to produce 
Dd- We can read from figure{T]that it seems the pole mass rric = 1.47GeV is more appropriate 
than the other two choices. The results are summarized in TableHl 

Now we turn to the analysis of 0"'' channel from scalar sum rule of Eq.f ll2p . The observed 
spectrum of 0'*' channel is reversed from the estimate from naive quark model. The numerical 
results following scalar sum rules of Eq.f ll2p are shown in figureEl We can see that under 
the threshold determined from Dg, a mass gap M^^ — M^^ ~ 15MeV between and Dg 
can be realized which is also lower than the experimental one ~ 35MeV. However, it indeed 
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TABLE I: Mass of pseudoscalar and Dg read from shaded areas marked by short bars in figure. [J 
for different threshold sq and rric at scale /U = IGeV. 



so(GeV2) mc(GeV) Dd{MeY) -Ds(MeV) M^^ - M^^^ (MeV) 
7.0 1.43 1812 1846 34 
6.8 1.47 1856 1888 32 
5.8 1.51 1873 1911 38 



TABLE II: Masses of scalar Dd and Dg read from figure. [2] for different threshold sq and nric at scale 
/i = IGeV. The first two values are read from the shaded area marked by short bars while the third 
is the central value between the two short lines in the third graph. 

so(GeV2) mc(GeV) Dd{Me\) -Ds(MeV) Mc^-MD,(MeV) 

7.9 1.43 2356 2340 16 

7.1 1.47 2341 2328 13 

6.5 1.51 2354 2340 14 



gives a correct splitting trend which agrees with experiment. The results in the O"*" channel are 
summarized in Table HTl 

It is instructive to study the scale dependence of our results since physical quantities are 
scale-independent thus it will supply a natural check on our results. Therefore the theoretical 
splitting should be unchanged when calculated with another scale. To this end we evolve the 
related parameters according to Eq. ffT^ to a higher scale /i = 1.3GeV which is still lower than 
charm mass. The results for 0~ and 0"*" are shown in figure. [3] and figure. H] respectively. We 
can see that when the scale increases, the results of 0~ channel still keep a well behavior as 
/i = IGeV. We can read from figure. [3] the splitting M/)^ — Mo^ ~ 35MeV which agrees well 
compared with = IGeV. But the situation is not so good in the 0"^ channel when we scale up. 
It is obvious from figure. H] that at /i = 1.3GeV the calculated mass is monotonically decreasing 
within the window satisfying the criterion 1 in 0^ channel. But fortunately the calculated mass 
of Dct is always larger than Dg within the selected Borel windows. For example, if we take the 
central values of in figure. H] we find the mass gap M/)^ — M^^ ~ 15MeV which is also consistent 
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1.90 - 
1.88 - 



0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2.0 

M''(GeV') 



FIG. 3: Mass of ^'^(dashed line) and Dg meson (solid line) of channel from pseudoscalar sum rule 
of Eq. (jl2p as function of Borel momentum Af^ at scale ^ = l.SGeV. 



with the sphtting obtained at = IGeV. Therefore the sphtting in both channels are invariant 
which shows a correct scale invariance. Combined the results at /i = IGeV we conjecture the 
reason why it is difficult to develop an extremum value at large charm mass in 0"'' channel 
maybe, as pointed out in23[, is that a large charm mass will induce large error. Our results 
imply that it is more appropriate to take a lower pole mass for charm. In fact as the scale 
increases, we approach to the asymptotic free side further thus the non-perturbative effects will 
have reduced impact. We can see obviously from figure JH that at the high energy side in the 
0"*" channel the mass gap decreases. 
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2.54 p 
2.52 - 




2.28 I I 

1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2.0 2.1 2.2 2.3 2.4 25 2.6 

M''(GeV') 



FIG. 4: Mass of Drf(dashed line) and Dg meson (solid line) from scalar sum rule of Eq. ()12l) as function 
of Borel momentum Af^ at fi = l.SGeV. It is obvious with the charm mass increasing, it is difficult 
to develop an reasonable extremum value in the allowed Borel window. 

The different importance of mass effects in realizing the splitting in 0~ and O"*" channel 
of D meson is not surprising. One should notice the "force" induced by the QCD vacuum 
or equivalently, the non-perturbative effects is parity-dependent which is well indicated by 
the contribution of the dominant condensates mc(gg) and mc{qcrGq) as well as the parts in 
(asG^/yr) introduced by their mixing. In 0~ channel the entire effect of these two terms give 
positive contributions to the correlation function thus an attractive force is induced by the QCD 
vacuum. While in O"*" channel these two terms supply negative contributions to the correlation 
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TABLE III: Mass of pseudoscalar and Dg read from shaded area marked by short bars in figure. [3] 
for different threshold so and rric at scale ^ = l.SGeV. 



so(GeV2) 


mc(GeV) 




I?s(MeV) 


Mb, -M^,(MeV) 


7.6 


1.43 


1807 


1845 


38 


7.2 


1.47 


1851 


1886 


35 


6.6 


1.51 


1885 


1922 


37 



function, then a repulsive force is induced by the QCD vacuum. Since —{dd) > —{ss) we could 
expect Md^ > Mo^ in 0~ and Mo^ < in O"*" channel. Furthermore the "force" is scale 
dependent which implies that the larger the scale is, the farther we leave from the confinement 
sector, therefore the importance of non-perturbative effects will be discounted compared with 
a lower scale. In the 0"*" channel the effect of the quark condensate (qq) overpowers other mass 
effects, so we find M^,^ < M^^. On the contrary if we set (dd) = (ss) then the mass difference 
in 0+ channel is produced by m^-dependent terms only; thus it is expected there will be mass- 
flipping. The mass curves are shown in figure. [51 It is obvious that the mass gap of the two 
scalars are very sensitive to the ratio k = {ss)/{uu). However, in 0~ channel the sign of the 
splitting remains unchanged, the mass gap of the two states is not sensitive to this ratio. 

In fact we can categorize corrections into two parts: one is parity-dependent and mainly 
proportional to light quark masses, {qq) and {qaGq), another is parity-independent such as 
'mq{qq), mq{qaGq) and (gg)^. It seems the former overpower the latter because the latter are 
doubly suppressed by the and the {qq)., however, their magnitudes for SU(3)-breaking are 
comparable. The former change their signs when we alter from 0"'' to 0~ channel and vice versa 
while the latter do not. We can learn from pseudoscalar sum rule in Eq.( lT2l) that in the 0~ 
channel these two parts provide a consistent response to the splitting since there is flipping in 
their signs. For instance the mass gap from quark condensates is: 

mc{{ss) - {dd)) > 0, 
md{dd) — ms{ss) > 0, 

and from the other condensates all contribute positive differences, all these positive difference 
are inclined to broaden the mass gap between Dg and in 0~ channel. However, in the O"*" 
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FIG. 5: Mass curves of (dashed line) and (solid line) vs. Borel momentum in 0"'' channel at 
= IGeV, where we set (dd) = (ss). It is clear the splitting of the two states reverse. 

channel all the non-perturbative corrections keep the same positive parity, and consequently a 
competitive role to the splitting of and Dg developed: 

mc{{dd) - {ss)) < 0, 
md{dd) — ms(ss) > 0, 

so does the (qcrGq) . This means there are some compensations to the splitting from the parity- 
independent terms which prefer to weaken the splitting between and Dg in O"*" channel. 
Therefore the mass gap is broader in the 0~ channel than in 0"^ of D meson. It is worthy to 



mention that this phenomenon is partly noticed in 23[ . In fact these results can be generalized 
to the 0~ and 0^ channels of pure-light mesons, but the mass difference induced by various 
condensates is greatly suppressed by small light quark mass thus it is expected the splitting is 
tiny by this mechanism, so to realize a realistic splitting in QCD sum rules based on a naive 
quark model should take into account instanton effects |8|, |9|. 

As operator mixing changes the coefficient of two-gluonic condensates significantly, we can 
see the mg-dependent parts in Cq^ are parity-dependent and have a complicated form but also 




FIG. 6: Mass curves of 0~(left) and 0"^ (right) channels vs. Borel momentum at fi = IGeV, where 
we have turned off the mg-dependent corrections in Cq2 in sum rules. The dashed and solid line in 
both graphs denote and Ds respectively. 

an obvious scale-dependence. It is instructive to investigate the impact on the calculated mass 
of both channels. For this purpose we turn off m^-dependent parts in Cq^ in sum rules. The 
mass curves are shown in figure|6l It is obvious that in this case in both channels Md^ is 
larger than M^^. We find in 0~ channel the mass gap is broader than there is mg-dependent 
corrections in gluonic condensate ~ lOMeV so our results show this correction is negative in 
the mass gap for the 0~ channel. In 0+ channel the situation is just opposite to the case when 
mg-dependent corrections turn on, thus the results show this contribution is positive to realize 
a reasonable splitting in 0''" channel. 

A natural idea is to generalize these arguments to the 1~ and 1+ channel of D. Unfortunately 
it does not work which can be well understood from heavy quark effective theory js^]. When 
there is no orbital excitation of light content, the c-quark with spin Sc = 1/2 and the light 
degrees of freedom with spin si = 1/2 forming a multiplet of hadrons with spin: 

J = ^®^ = 0©1, 

thus it is clear there is a unique = 1~ multiplet. While if there is orbital excitation of light 
degrees of freedom, the spin of light content will be: 

Si = I ± - = - or -, 
' 2 2 2' 
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combining with the spin of c-quark Sc = 1/2: 

so there will be two 1^ multiplets, the D* and the Di states experimentally. These two states 
are very close, thus the single resonance approximation in QCD sum rules is not viable. If 
we let the charm mass go to infinity these two 1"*" states can be separated in the formalism of 
the heavy quark effective theory 381]. However, the l/rric corrections which are the same order 
of SU(3)-breaking effects make these two states mix again. So we still cannot get a simple 
correspondence between 1~ and 1"*". We can resort to the experimental data directly 36|] . The 
mass splitting between Dg and in 1~ is about lOOMeV, while it is only about 40MeV in the 
lower l"''(which corresponds to 1~ in the heavy quark mass limit). A similar effect still appears. 

Finally let us briefly mention the B case although there is not enough experimental evidence. 
Since the b-quark mass is so large the SU(3)-breaking effects are smeared in the formalism 
Eq.f ll2p . We can hope the similar effects will be recovered in the formalism of the heavy quark 
effective theory. Certainly, SU(3)-breaking effects also appear in l/rrih corrections which might 
not be small and could cause some differences between D mesons and B mesons. 



IV. CONCLUSIONS 

In this work, based on the pseudoscalar and scalar sum rule from cq structure we inves- 
tigate the SU(3)-breaking effects enhanced by the large charm mass on the splitting of the 
pseudoscalar and scalar D multiplet. Since the quark condensates {qq) and (qaGq) are greatly 
enhanced by the heavy quark mass, they play more important roles in the SU(3)-breaking in 
the heavy-light mesons. The sign of the (qq) and (qcrGq) contributions is different in 0~ and 
0"^ channels thus resulting a " parity- dependent " force which can explain the relatively-lower 
mass of i5s(2317). Furthermore, these parity-dependent corrections broaden the mass gap in 
the 0~ channel while weakening it in 0"*" channel. The results show that the 0{mq) corrections 
in two-gluonic condensates introduced by operator mixing are noticeable in 0"'' channel. Such a 
" parity-dependent " force should exist generally in the heavy-light mesons, but its magnitude 
should be dependent on the specific system. This force is also energy-dependent, so its effect 
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should be weaker in high excited states. We also analyze the cases 1 and l"*". The experimental 
data hints such a force also exists but is weaker than that in 0~ and 0"*". 

V. ACKNOWLEDGEMENTS 

This work is supported by NNSFC under Projects No. 10775117, 10847148. The author 
H.Y. Jin thanks Hai-Yang Cheng for very useful discussion. 

Appendix A: derivation of Eg. (171) 

It is a little effort to work out Cg^ in scalar current expansion from a vector current expan- 
sion. To this end, let us consider vector current two-point function: 



n^.(9') = I j rf^xe*^^(0|r{g(x)7Mc(x),c(0)7.g(0)}|0) 

= {-9t.uq^ + q,qu)^^{q^) + g^g.nS(g2), (Ai) 

To single out the scalar part we contract Eq. llAll) with q^q'^: 

qy^HM) = ?'nS(g2), (A2) 

Then it is convenient to consider the following two-point function based on the four- divergence 
of vector current: 



= i j rf^xe^'''^(0||g(x)^c(x) + g(x)^c(x)} x |c(0)^g(0) + c(0)^g(0)}|0) 

= i(me-mg)2 [ rfVe*'?^'(0|g(a;)c(x)c(0)g(0)|0), (A3) 



where the Dirac equations: 

~^ip{x) = —imip{x) 

and 

- "tr 

ijj{x) = imip{x). 
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have been used. In this way we can associate OPE of scalar current expansion with the vector 
one: 

I [ tiW^^(0|g(x)c(x)c(0)g(0)|0) = , " (A4) 
J [nic - mqY 

Combined with Eq.( ]A2l) it is straightforward to obtain: 



i / rfV(0|g(x)c(x)c(0)g(0)|0) = (A5) 



{rric - mqY' 



And C^s has been worked out 39| in heavy-hght vector current expansion, here we write it in 
a standard form: 



where 



A = :--:(-30W^ -SW- 3), 



W 

2880^m6 

,2 



B = --{-lOW^ + m'^ + 31^ + 2), 



— 



Combined with Eq.f lASI) we can obtain C^a in a scalar heavy-light expansion: 



2 



= -^^^(-10^' + 4iy + 3iy + 2). (A7) 



This is Eq.([7]) where we have suppressed the factor 
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